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Abstract. A method for testing stability of relative equilibria in Hamiltonian 
systems of the form "kinetic + potential energy" is presented. This method 
extends previously existing techniques to the case of non-free group actions 
and singular momentum values. A normal form for the symplectic matrix at 
a relative equilibrium is also obtained. 
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1. Introduction 

A relative equilibrium in a dynamical system with symmetry is a point in phase 
space for which its dynamical evolution is contained in a group orbit. The study 
of relative equilibria in symmetric Hamiltonian systems has been around for a long 
time, with its origins in the field of analytical dynamics and more recently using the 
modern symplectic and Poisson geometric setup. Relative equilibria are important 
since they are the analogues to equilibrium states in systems with symmetry, for- 
malized with the action of a Lie group on the phase space. In physical applications, 
the only observable equilibrium states are those which are stable under small pertur- 
bations. Similarly, in the symmetric context, the only observable relative equilibria 
are those which are stable in some adequate sense. Based in Noether's Theorem, ge- 
ometrized in the property of the invariance of the level sets of the momentum map, 
the notion of stability generally adopted is that of G^- stability introduced in [15j . 
and that is closely related to the Lyapunov stability of the induced Hamiltonian 
flow on the reduced phase space. 

In the field of analytical dynamics, the classical Routh Theorem gives conditions 
on the stability of steady motions which keep stationary the value of a first integral 
of a dynamical system for fixed values of the others (see for instance [18] and the 
treatments based on the Routh Theorem in [4j[T4]). Relative equilibria are seen 
in this context as steady motions for systems having cyclic coordinates due to the 
existence of a symmetry group, for which the components of the momentum map 
together with the energy provide a set of first integrals. In the last decades, the 
implementation of these principles within the field of Geometric Mechanics has been 
studied. This has produced methods (like the Energy-Momentum Method p~5] and 
the Energy-Casimir Method pQ, see also [9] for an overview) to test the stability 
of relative equilibria in Hamiltonian systems for arbitrary symmetry groups and 
momentum values. These methods exploit Noether's Theorem and the symplectic 
and Poisson geometry of the phase space. In the case that the relative equilibrium 
under study lies in a regular value /i of the momentum map, the Energy-Momentum 
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Method of [15] provides a technique to test its stability modulo the action of G^, 
the stabilizer of the momentum value [i under the coadjoint representation of G. 
This was generalized in [5] and to also cover the case when the momentum 
value is singular (this happens when the symmetry group does not act freely on 
phase space), assuming G^ is compact. Also, in [10| and [16] stability of relative 
equilibria satisfying several other hypotheses is investigated. 

A very important particular kind of Hamiltonian system is the class of simple 
mechanical systems, paradigmatic of Classical Mechanics, since many Hamiltonian 
system of physical interest lie in this category or can be obtained from a simple me- 
chanical system by some suitable reduction process. These have as phase space the 
cotangent bundle of a Ricmannian manifold (called configuration space) equipped 
with its canonical symplectic form, and the Hamiltonian function is of the type 
"kinetic plus potential" energy, where the kinetic energy is given by the norm ob- 
tained from the Riemannian metric, and the potential energy is the pullback of a 
function defined on the base. Symmetry in this systems is implemented by the lift 
of an isometric action on the base that preserves the potential energy. This big 
amount of extra structure with respect to general Hamiltonian systems on arbi- 
trary symplectic or Poisson phase spaces implies that in simple mechanical systems 
everything is constructible from the knowledge of the configuration space, its Rie- 
mannian structure, the action of the symmetry group on it and the choice of an 
invariant potential energy. Therefore, it is reasonably to expect that the stability 
methods referred previously will particularize in a way that the involved computa- 
tional complexity will simplify considerably, in that one would work at the level of 
the configuration space, instead of on its twice dimensional cotangent bundle. 

In the case of regular relative equilibria, this refinement of the Energy-Momentum 
Method has been worked out in [19], and the obtained stability test for relative 
equilibria in simple mechanical systems is known as reduced Energy-Momentum 
Method. Its conditions for G M -stability are reduced from the level of phase space 
to the level of configuration space. This method has the highest degree of sophis- 
tication among the different stability tests available in the literature of symmetric 
Hamiltonian systems, and as part of it, it provides a block-diagonalization tech- 
nique that allows to express the linearization of the Hamiltonian vector field at 
a relative equilibrium in a way adapted both to the symmetry of the system and 
to the fibered structure of the phase space. This block-diagonalization yields also 
further simplifications in the stability analysis. 

Surprisingly, in the very frequent and important case of singular momentum 
values such a refinement for simple mechanical systems has not been studied in 
detail, and thus the application range of the reduced Energy-Momentum Method 
is severely limited. Indeed, the literature of applications of the theory of relative 
equilibria is full of examples in which singular relative equilibria of simple mechan- 
ical systems are studied with general geometric and Hamiltonian techniques which 
neglect their extra structure, in particular for the stability analysis. This paper 
provides a solution to this situation by obtaining a generalization to singular mo- 
mentum values of the reduced Energy-Momentum Method and its main features. 

In Section 2 we quickly review the theory of relative equilibria for general Hamil- 
tonian systems and simple mechanical systems, and we collect some of the standard 
results on their stability by geometric methods. Section 3 is a necessary technical 
interlude on the properties and geometry of a distinguished symplectic component 
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of the linear slice for a cotangent-lifted action, and most of our subsequent results 
will rely on this section. In Section 4 our main result, Theorem 14. 1\ is stated, pro- 
viding an extension of the reduced Energy-Momentum Method of 19J applicable 
to relative equilibria with singular momentum values. Section 5 applies this result 
to a classical example of a relative equilibrium with a singular momentum value in 
a well-known simple mechanical system consisting of an axisymmetric rigid body 
with a fixed point in an homogeneous gravity field. It is shown how the application 
of our method simplifies the stability analysis with respect to the application of 
the methods developed for general Hamiltonian systems. In Section 6 we extend 
the block-diagonalization result of [19] to the singular case, in Corollary 16.21 and 
Proposition 16.21 Finally, Section 7 puts in context our results with related work 
in the literature. In particular it is shown how the block-diagonal expression for 
the symplectic matrix of Proposition 16.21 particularizes in the regular case to the 
normal form obtained in [19j . and a comparison is also made between our results 
and the Lagrangian Block-Diagonalization method of Lewis [6] . 

2. Relative Equilibria and simple mechanical systems 

Let (V,lu) be a smooth finite dimensional symplectic manifold with symplectic 
form u and G a finite dimensional Lie group acting smoothly, properly and in a 
Hamiltonian fashion on (T 3 ,^) with Ad*-equivariant momentum map J : V — > g*. 
Given a G-invariant Hamiltonian function h £ C G (V), a point z £ V is called a 
relative equilibrium (for h) if its Hamiltonian evolution lies inside a group orbit. 
Equivariance of J and Noether's Theorem imply that the Hamiltonian evolution 
of z is described as the orbit of z by a one-parameter subgroup of G generated by 
a Lie algebra element £ which belongs to g M C g, where /i = J(z) and g^ is the 
Lie algebra of G M , the stabilizer of /i for the coadjoint representation of G. The 
element £ is called a velocity of the relative equilibrium. Using the usual notation 
for the infinitesimal action of g on V the condition for z to be a relative equilibrium 
is written as Xh(z) = £-p(z), where Xh is the Hamiltonian vector field associated 
to the function h. 

If the stabilizer G z of z is not discrete then there is a degeneracy in the choice 
of a velocity for a given relative equilibrium, since any representative of the class 
[£] £ g/g z produces the same orbit of z. In any case, by equivariance of J, the 
inclusion G z C G p holds. The quintuple (V, u>, G, J, h) will be called in short a 
symmetric Hamiltonian system. 

The following definition introduced in 15] is generally adopted as the correct 
notion of stability of relative equilibria in Hamiltonian systems, generalizing in the 
Hamiltonian context the concept of Lyapunov stability of fixed equilibria for flows 
of vector fields. 

Definition 2.1. A relative equilibrium z with momentum [i — J(z) is said to 
be G ^-stable if for every G ^-invariant neighbourhood U of G M • z there exists a 
neighbourhood O of z such that the Hamiltonian orbit of O lies in U. 

In [5J 111] a method for testing stability of relative equilibria with singular mo- 
mentum values in Hamiltonian systems is developed, generalizing the Energy- 
Momentum Method of [15] for relative equilibria with discrete stabilizers. We 
quote here the main result, due to its importance in the subsequent development 
of the paper. For that, given an element £ £ 9, define the augmented Hamiltonian 
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hi e C°°{V) as 

(2.1) h i (z)=h(z)-{3(z),0- 

It well-known that z G V is a relative equilibrium for the symmetric Hamiltonian 
system (V, u, G, J, h) with velocity £ if and only if z is a critical point of h^. Also 
recall that since the Hamiltonian G-action is proper, any stabilizer for this action, 
in particular G 2 , must be compact (see [2]). Since q z C g M we can choose a G z - 
equivariant splitting of as g p = q z © g^" . 

Theorem 2.1 ([5]. [11]). Lei (V,Lu,G,3,h) be a symmetric Hamiltonian system 
and z €z V a relative equilibrium with stabilizer G z and velocity £. Assume that 
J(z) = /i and that G^ is compact. Then £ 6 0^. Lei £ &e t/ie projection of £ 
onto some G z -invariant complement of q z in g M (always available by compactness 
of G z ). If & 2 z h^±\ v is definite for some (and hence any) complement V s to g M • z in 
kerT z J, then z is G ^-stable. 

In this theorem the ambiguity in the velocity introduced by the stabilizer of 
the relative equilibrium appears explicitly. In typical computations, testing this 
condition over all possible G z -invariant complements of q z in g M gives the sharpest 
stability results (see [11] and the example in Section 5) . There is an infinite number 
of choices for the space V s in Theorem 1 2. 1[ and any of them is called the (maximal) 
symplectic normal space at z, since it is a maximal symplectic subspace of the 
symplectic orthogonal to the group orbit at z. 

In this paper we will study a particular case of Hamiltonian systems of great 
interest in Classical Mechanics. This is the class of the so-called symmetric simple 
mechanical systems, which are symmetric Hamiltonian systems where V is T*Q, 
the cotangent bundle of a smooth, finite-dimensional Ricmannian manifold (Q, -C 
•, • 3>) equipped with its canonical symplectic form lu and G is a finite-dimensional 
Lie group acting properly and isometrically on Q and by cotangent lifts on T*Q. 

Following [20] , the Hamiltonian function is constructed in the following way: let 
V be a smooth G-invariant function on Q and call V — t*V £ C G (T*Q), where 
r : T*Q ~ * Q is the cotangent bundle projection. We will refer to both V and V 
as the potential energy. Let FL : TQ — > T*Q be the Legendre map associated to 
<C ■, ■ defined by the formula 

(2.2) (FL(v x ), w x ) =<C v x , w x Vv x ,w x eT x Q. 

The Legendre map is a G-equivariant vector bundle isomorphism covering the iden- 
tity on Q. With it, we can define the kinetic energy K £ C (V) as 

K(p x ) = \ < FL- 1 (p a ),FL- 1 ( fe ) », V Px £ T*Q. 
Finally, the Hamiltonian h is defined by 

(2.3) h = K + V. 

With respect to the canonical symplectic structure of T*Q the cotangent-lifted 
action of G is Hamiltonian, with equivariant momentum map defined by the ex- 
pression 

(2.4) (J(p*),0 = (p*,£q(z)>, V£e . 

The symplectic manifold T*Q is the phase space of the Hamiltonian system, while 
the base Q is usually called configuration space. Accordingly, for any point p x in 



STABILITY IN SIMPLE MECHANICAL SYSTEMS 



5 



T*Q, the projection x = r(p x ) is called the configuration point (or base point) of 

Px- 

A key feature of simple mechanical systems is that both their geometric and dy- 
namical properties are entirely constructed using the knowledge of the Riemannian 
manifold (<3,<§; •, • 3>), the isometric G-action on it and the choice of a potential 
energy V. Thus, one could reasonably expect that the implementation on this class 
of Hamiltonian systems of the stability test given in Theorem 12.11 should simplify 
accordingly and yield easier computations at the level of Q and G instead of the 
bigger space V = T*Q. The obtention of such a refinement to simple mechanical 
systems of this stability test is the main result of this paper, Theorem 14. 11 

Remark 2.1. There are several treatments of this problem in the literature. Simo 
et al. develop in [T9] a sophisticated test, that particularizes Theorem 12.11 called 
the reduced Energy-Momentum Method. This test gives sufficient conditions for 
the stability of relative equilibria in simple mechanical systems provided the config- 
uration point of the cotangent relative equilibrium has a discrete stabilizer (which 
is the same as to require that its momentum value is regular) . The main advantage 
of this is that it is constructed specifically for this class of systems, and this fact 
reflects in less computational difficulties than the application of the main method, 
Theorem 12.11 designed for general Hamiltonian systems. In [6 a Lagrangian ana- 
logue of the results of [19] is obtained, being valid also for relative equilibria of 
a larger class of mechanical systems. Here, based in Theorem 12.11 we produce a 
method for testing stability of relative equilibria in simple mechanical systems that 
could be seen as a generalization of the reduced Energy-Momentum Method to the 
singular case, i.e. without requiring discrete stabilizers of configuration points or 
regular momentum values. 

3. A COTANGENT-BUNDLE ADAPTED SPLITTING OF THE SYMPLECTIC NORMAL 

SPACE 

In this section we describe a realization V s C T Vsc (T*Q) of the symplectic normal 
space at a relative equilibrium p x of a simple mechanical system, as well as a 
convenient cotangent-bundle adapted splitting of V s which will be extremely useful 
for the remaining constructions. Most of the results of this section are merely 
expository, and a complete description including proofs and the obtention of the 
symplectic normal space at points p x of general form, not only relative equilibria, 
can be found in [T?] . 

One of the geometric objects that will be extensively used is the locked inertia 
tensor I, a family of bilinear positive semi-definite symmetric forms on q defined 

by 

(3.1) I(a?)K,»7) =<& (*), Ve.rjGfl, xtQ. 

Note that at each point x, the kernel of I(x) is precisely g x , the Lie algebra of the 
stabilizer of x. Therefore i(x) is a well defined inner product on g only at points of 
Q where the action is locally free. The locked inertia tensor satisfies the following 
invariance and infinitesimal invariance properties (see [9]): 

(3.2) \{g-x){KA g ^kA g i 1 )=l{x){^ 1 ) 



(3.3) (DI ■ Xq(xM,v) +I(x){adxt,v) +I(x)^,ad\v) = 0, 
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for every g £ G, x £ Q and £, 77, A € g. Note that G Px C G x by the equivariance of 
t. Let a; € Q be the base point of an element p x E T*Q and denote by H = G x its 
stabilizer. Using the G-invariant Riemannian metric on Q we can form the splitting 

(3.4) T X Q = g ■ x © S 

where S = (g ■ x) , and it is usually called a linear slice (for the G-action at x). 
Hence, S is the space of directions orthogonally complementary to the group orbit. 
This is obviously a //-invariant splitting for the induced linear //-action on T X Q. 
Next we choose a G Px -invariant splitting of the Lie algebra g of the form 

(3.5) = f)ffit. 

This is always possible by the properness of the G-action on Q, which implies that 
G Px is compact. A concrete way of choosing (|3.5p will be introduced in (|3.14j) . For 
any element (eg we write in a unique way £ = + £ r , relative to this splitting. 
The space r collects the elements of g that generate nontrivial orbits of x. 

Noting that r ~ g ■ x by the isomorphism £ 1— > £,q(x), we can compose this 
identification with (|3.4|) and dualize, to get 

(3.6) T x Q~t©S and T;Q~r*©S*. 

Associated to the Riemannian structure on Q, there is an Ehresmann connection 
on T*Q, for which the connection map at p x , K : T Px (T*Q) — > T*Q is defined as 

(3-7) K(X) = 



f= c(t) VXET Px (T*Q), 



the covariant differential of c(t) along c(t) = r(c(t)) relative to the Levi-Civita 
connection V. Here c(t) is any local curve c : (— e,e) — > T*Q projecting to c(t) 

and satisfying i) c(0) = p x , and zz) t Q c(t) = A. This connection map K at 

Pa is a G Px -equivariant linear map, which combined with the differential at p x 
of the cotangent bundle projection r yields a G P:c -equivariant linear isomorphism 
* : T Px (T*Q) -> T a Q © T^Q = T*(T X Q) defined by 

*(A) = (T Px r(A),A(X)) VAeT Px (T*Q) 

(see [17]). We call vectors at p x lying in the kernel of T Px r vertical, since they are 
tangent to the cotangent fiber through p x . Analogously, those elements of T Px (T*Q) 
which are in the kernel of K are called horizontal, and are identified through T Px r 
with vectors tangent to Q at x. We now compose the above isomorphism ^ with 
the two dual isomorphisms (|3.6|) to get a new one 

(3.8) I:T Px (T*Q)^(x®S)®(x*®S*) 

which can be explicitly expressed as I{X) = [rj,a;v, a), for the unique r),a,v,a 
satisfying 

t p* t ( x ) = a + VQ{x) 

K(X) = a + ¥L((^ 1 (p)) Q (x)y 

Here lo denotes the restriction of I to r, according to (|3.5p . Note that now So 
becomes a well-defined inner product in r and thus also a linear isomorphism Iq : r — > 
r* ~ t)°. We can therefore work in the image of /, which we call /-representation, 
instead of on T Pj .(T*Q), and that is what we will do in the rest of the paper. Note 
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that in this identification the space of vertical and horizontal vectors is expressed, 
respectively, as 

(0,0; v, a) VaeSVet*, and (77, a; 0,0) Va e S, 77 e t. 

The isomorphism / is G Px -equivariant with respect to the linear action on the target 
space given by 

(3.9) g- (ri,a;v,a) = (Ad 5 77, g ■ a, Ad*_i v, g ■ a), 

where g ■ a and g ■ a denote, respectively, the restriction to G Px of the linear repre- 
sentation of H on S and its contragrcdient representation on S*. This action is well 
defined since G Px C H and S, r, and their duals are //-invariant by construction. 

In order to obtain a convenient characterization of the symplectic normal space 
V s C T Px (T*Q) at a relative equilibrium in the /-representation, and also for future 
reference, we quote some technical results introduced in [17) . There it is proved 
that it is possible to extend vectors v £ S to local vector fields v defined in a 
neighbourhood of x, in a way adapted to the G-action, and such that the family of 
vector fields Xq,v, for any A £ 2 and v £ S spans T X >Q at every x' near x. 

We will sketch here the obtention of the local field v. Recall that by Palais' 
Tube Theorem [12] we can construct an invariant tubular neighbourhood of an 
orbit G ■ x c Q as follows: Let H = G x act on G x S as h ■ (g, s) — (ghr 1 , h ■ s). 
Let G x h S be the quotient space for this action. Then there is a H- invariant open 
ball U C S centered at the origin such that the map a : G Xh U — > Q defined by 

(3.10) cr([g,s}) = g ■ exp x s 

is a diffeomorphism onto a G-invariant neighbourhood O of G-x. Here exp^ : T X Q — > 
Q denotes the exponential map associated to <C •, • This diffeomorphism is G- 
equivariant with respect to the given action on Q and the G-action on G Xjj S 
defined by g' ■ [g, s] = [g'g, s]. 

Now choose any inner product on g such that the splitting (|3.5[) is an orthogonal 
direct sum. Extend this inner product by right translations to a //-invariant Ric- 
mannian metric on G. Then we can interpret r as a linear slice at the identity for 
the free //-action on G given by (h,g) 1— > gh^ 1 . It follows that if we call exp e the 
exponential map for this metric on G, there is a small //-invariant neighbourhood 
of e in G such that every g belonging to it can be written as g = exp e £ r /i -1 for 
unique elements £ r £ r and h £ H. Using the Tube Theorem every element in a 
small neighbourhood of x in Q (not G-invariant in general) can be expressed as 
x' = er([exp e £ r /i -1 , s]) for unique elements £ r £ r, h £ H and s £ U. It is easy to 
prove that for any x' = cr([exp e ^h^ 1 , s}) near x, the formula 

F*(x') = a([ex Pe f fe -1 , a + • «]) 

defines a flow Z 7 ^ for any v £ U C S. The associated local field u is then defined as 

(3.11) ^x>) = ^ o Kix>) 

for every x' near a;. In the following theorem we collect the most important prop- 
erties of this family of local vector fields. For the proof, see [17] . 



Theorem 3.1. Letv,w£ S and 77, A, £ g. Th 

(ii) « V^^iQ^), W>=-\ (DI • «,) (^, 
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(iii) « V Cq v(x) , Xq (x) »= \ (DI • v) A) 

(iv) « V-£ Q (x), Xq(x) »= | (DI • v) (£, A) 

(v) < V-£ Q (x), w >=< V ?Q u(a;), w > + < £"> • v, w >. 

Here, and ^ denote the projections of elements of g onto r and [j according 

to (1331) . 

Notation: We will introduce for any i> € S a linear map G(v) : r — > S* defined as 
(3.12) < C(«)(D, ^ >s=< V ?Q u(a;), u- >, 

where <C •, • is the restriction of the metric on Q to S C T X Q. Note that G is 
not linear in v since it depends on the concrete extension v. We will also employ the 
following notation: if W is a linear subspace of the linear space V and i : W ^ V 
its inclusion, we will write Pw : V* — > VF* for its dual projection. 

We fix from now on a point of the form p x = FL(£q(x)) with G x — H . It follows 
from (|2.4p and the definition of the locked inertia tensor that p x has momentum 
/i = J(p x ) = I(.t)(£) = Io(£ c ). The reason for this choice of p x will be clear in the 
following section, where it is explained why every relative equilibrium of a simple 
mechanical system must be precisely of this form. 

Remark. It is a well-known property of points of the form p x = FL(£q(x)) that 
G Pl , — H n G M . This follows immediately from the relation G Px = H Px and iden- 
tifying p x = (0, /i) <G S* © r* using the //-isomorphism (|3.6j) . For general points of 
T*Q one has only an inclusion G Px C H n G p . 

We now make a concrete choice for the complement r in (|3. 5|) , as well as for other 
relevant subspaces of g. Start by choosing a G Px -invariant complement p to g Pa . in 
M , i-e. 

(3-13) Al =0 P . OP- 

Next, define a G Px -invariant complement t to t) © p in g in such a way that defining 
t = p © t we have that p and t are orthogonal with respect to the restricted locked 
inertia tensor lo . We can then write 

(3.14) g = f)0t=f)0p©t, 

and hence we have constructed the splitting (|3.5p . 
Let us define the following subspace of g 

(3.15) ^ = {A e t : P„ [ad>] = 0}. 

This space will play an important role in our characterization of V s , and it can be 
proved (see |17j ) that it is isomorphic to the symplectic normal space at /i for the 
restriction to H of the coadjoint action of G on O^, the coadjoint orbit containing /i. 
Note also that by using {331) w e can write q M = {A e t : P|, [(DI • A Q (.x))(<f ) = 0]}. 

As a particular case of Theorem 6.1 in [17] the space kerT Px J consists in the 
/-representation in the elements (77, a; v, a) € (r S) © (r* © S*) satisfying 

(u, A r ) - \ {(DI ■ i Q [x)) (X\ rf) - (DI • a) (A r , £) - (ad^/i, ry)} + (ad* h/ i, ry) = 
for every A r S r and A 1 ' S f). Also, in the /-representation 

0m-P*= |(A,0;ip t [(DI^ Q ( a ;))(A)] ) -ip s [(DI.(-))(A,O]) : VAepj. 
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From the above two expressions it is easy to obtain that the symplectic normal 

space V s , a complement to g M • p x in ker T Vsc J , can be chosen to be 

(3.16) 

V s = {(A,a;iP r [(DI-^(x))(A)+ad^-(DI-a)(e)],7 

- iPs [(DI-(-))(r,A)] + «C(a)(r),-»s) : VAeqf, aeS, 7 eS*}, 

with q M defined in ()3.15j) and S in (|3.4j) . The symplectic normal space V a is G Px - 
invariant by construction with respect to the action (|3.9|) in its ambient space (see 

4. Stability of singular relative equilibria in simple mechanical 

SYSTEMS 

In the following we will be in the setup of Section [5] and fix a simple mechanical 
system h = K + V as in (|2.3[) . In this framework, once an element £ G g is chosen, 
we can separate the augmented Hamiltonian (|2.1jl into a kinetic and a potential 
part as 

h(- = + V(: , where 

K £ (p x ) = l\\ Px -X £ (x)\\ 2 and V s (x) = V(x) - ±I(x)(£, £)• 

As for the potential energy we have used the notation V£ = r*Vf . The one- form x 
is defined by 

(4.1) x£(x) = FL(£ Q (x)). 

The functions and are called the augmented kinetic energy and augmented po- 
tential energy respectively. Recall now (see for instance Theorem 4.1.2 and Propo- 
sition 4.2.1 in |9J) that with the introduction of these two auxiliary functions we 
have the following characterization of relative equilibria: 

Proposition 4.1. Let p x G T*Q. The following are equivalent: 

(i) p x is a relative equilibrium of (|2.3[) with momentum /i and velocity £ G g^ 

(ii) p x is a critical point of the augmented Hamiltonian 

(iii) p x is simultaneously a critical point of and 

(iv) p x = FL(£q(x)) and x is a critical point ofV^. 

Note that (iv) restricts the form of phase space points candidates to be relative 
equilibria of (12. 3|) to be of the form p x = FL(£q(x)). Thus is the reason for studying 
in detail in the previous section the symplectic normal space only at this class of 
points. 

Let p x = FL(£q(x)) £ T*Q be a relative equilibrium for the simple mechanical 
system (|2 .3[) with momentum J(p x ) — /i. We call H = G x and we choose a 
(G Px = H n G M )-invariant inner product on g relative to which we construct the 
splittings (|3.13p and (|3.14p . Note that by hypothesis x is a critical point of 
for any £' G g^ such that [£' — £] = G g^/g P:c - In other words, any such £' is a 
velocity for the relative equilibrium p x . In particular, this happens for £ G p, the 
projection of £ onto p according to (|3.13|) . 

Let 5p G T P:c (T*Q) be a tangent vector at p x . We will write its horizontal 
and vertical components as 5p H = T Pa: T(8p) G T X Q and 8p v — K(Sp) G T*Q 
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respectively. It is clear that if I(Sp) = (A, b; u, 7) is the /-representation of Sp one 
has that 

Sp H = X Q (x) + b 

SP V = FL((I - 1 (i/)) ( a; ))+ 7 . 

Also, using (|3.6|) we can express the horizontal and vertical variations 5p H and Sp v 
as elements of r S and r* S* respectively like 

S P H = (A, 6) 

We will use both notations indistinctly. 

Finally, for a curve c(t) 6 Q with c(0) = a; we will write Hor Pa .(c(i)) for its 
horizontal lift to T*Q at the point p x with respect to the Levi-Civita connection. 
Equivalently, Hor Px (c(<)) is the parallel translation of p x along the curve c(t). 

Lemma 4.1. Let p x = WL{^q{x)) £ T*Q be a relative equilibrium for the simple 
mechanical system (|2.3|) with momentum /i and velocity £ £E Let £ the orthog- 
onal projection of £ onto p according to (|3.13[) . Then, for any 6pi, Sp 2 G T Pm (T*Q) 

(i) dl\%I(S Pl ,S P2 ) = d£Vt±(8p?,8p?) 

(ii) d 2 Px K^ (S Pll S P2 ) =< ^ - IT*X^ ■ 5p^) v ,5pX - (T xX * X ■ 5 P %) V ». 

Proof, (i) follows immediately by the definition of Sp H and noting that V^j. = 
t*V^±. To prove (ii) we will consider horizontal and vertical vectors separately. Let 
I(5pi) = (0;5pY) for i = 1,2. Then 

dlK^(5 Pl ,5p 2 ) = lA\ s = o A\ t ^Jp x + t5pY + sSpX-x^(x)f 

= llL= i| t =oii^+ s ^)ii 2 

= <t:5pY,s P Y -» . 

If /(<5pi) = (0;5p] / ) and I(Sp 2 ) — (Sp 2 ;0), let c| piJ -(x) be any smooth curve satis- 
fying c° p „ (x) = x and = c^ pf [x) = 6p% . Then 

d^x(<Jp 1)( Jp 2 ) = |^L = || t = II Ror px+tSp v(c° p¥ (x))-rf X (c° p „(x)) || 2 
= -i\ t = o^^s P ^X i± {x),p x +t5pY ~x e (x) » • 

But since p x is a relative equilibrium with velocity £, then p x = \^ ( x ) and then 

dl x K c ±(6 Pl ,Sp 2 ) = - < V^x^x),^ > ■ 

Finally, consider two variations of the form I(5pi) = (Sp^; 0) and L(Sp 2 ) = (<5p^; 0). 
Then 

d 2 px K i ±(Sp 1 ,Sp 2 ) = 

= hll = osL = II Hor Ho tpi ( c lf (,)) (Fjs (*)))- x s± (4 pf (*))) II 2 

= -lit = << (4 pf (x)), Hor^ (c* pf (*)) - (c* pf (zj) » 

=< V Sp ?X $± {x),V Sp? x $± (x) > ■ 

Where if Sp 2 — (A, 6), then 8p^} is the local vector field b + Aq, and — denotes 
its flow. Recalling now the definition (|3.7p of the operator if, for any v £ T X Q and 
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local extension v one has 



V„ X « (x)=\7 vX i (x)=K(T xX t -v) = (T xX t ■ v) 



and the result is proved. 



□ 



Since every horizontal variation 5p H can be written as the sum of two contri- 
butions one coming from r and the other from S, that is Sp H = Cq( x ) + b, for 
C S t, b G S, we can consider these two contributions separately and thus obtain 
concrete expressions for (T xX ^ ■ Sp H ) v . 

Lemma 4.2. Let ( e t and b £ S, and identify T*Q with r* © S* by the isomor- 
phism (|3.6p . Then, 



n (T xX £X <Q(x)) V = 

(±P r [(DI.^(x))(C) + ad^ + 2(DI.C Q ( a :))(^)],-iPs [(DI • (•)) (£\ C)] ) • 
(ii) (T xX ^ ■ b) v = (|P r [(DI • b) (^)] ,« C(b)(e\ ■ »s). 



Proof. The proof is just an immediate consequence of Theorem 13.11 in particular 
of items (i), (ii), (iv), and (v). We will just prove (ii). Recalling that V is a metric 
connection, then 



Similarly, if w G S, then by item (v) of the above referred theorem we have 



This yields result (ii) The proof of (i) is identical, including some manipulations 
using the infinitesimal equivariance property of the locked inertia tensor given in 



We apply now the results obtained so far in order to produce a singular version 
of the reduced Energy-Momentum Method of [19]. Consider V = T*Q with its 
canonical symplectic form in the statement of Theorem l2.il where (Q,<C •, • 3>) is 
a Riemannian manifold on which the Lie group G acts by isometrics. The Hamil- 
tonian action on T*Q is the cotangent lift of the action on Q, and the Hamiltonian 
is given by (|2 . 3[) defining a simple mechanical system. Let £ be an element of the 
Lie algebra g of G. Fix p x = FL(£q(x)) with momentum [i. Assume that 

(i) x is a critical point of , and 

(ii) G^ is compact. 

Then p x is a relative equilibrium for our simple mechanical system with velocity 
£ satisfying £ £ and Theorem 12.11 can be applied to study its G^-stability. 
That is, we need to determine when d^/i^i| y is definite. For that, we will use the 
characterization of the symplectic normal space V s , given in (|3.16|) which establishes 
a linear isomorphism k : q M ® S S* -> V s C T V T*Q ~ (rffi S) (r* S*) explicitly 
expressed as 




((T xX ^ -b) v ,X Q (x)) =«V F £^),A Q (z)»= i(DI-6)(^,A). 



((T xX « • b) v ,w) =« Vf$(x),w »=« C(b)(£ x ),w » s . 



(El 



□ 



(4.2) 



k (A, a, 7) = (A, a 



±P r [(DI • ^(z))(A) + adl M - (DI • a)(£ x )] , 
7 - iPs [(DI • (■))(^ ± , A)] + « C(a)(^), • »s) • 
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This map is G Px -equivariant with respect to the action on q M S © S* given by 

9 ■ (A, a, 7) = (Ad s A, g ■ a, g ■ 7), 
and the action $379j) on (r S) (r* ffi S*). 

Definition 4.1. At a relative equilibrium p x — ¥L(£q(x)), with momentum /i 
velocity £ € g M , define the linear subspace £ of T X Q isomorphic to q' 1 S, as 

(4.3) V = {\ Q (x) + a € T X Q : A e q^, a e S} , 
where q M is defined in p.!5|) . 

for any fwo vectors Vi,i)2 € T X Q, define the correction term as the symmetric 
bilinear form on T X Q defined by 

(4.4) corr^Xwi,^) = (P t [(DI ■ Jo 1 (P« I(DI ■ «a) «)])>■ 
We can now state the main result of this section. 

Theorem 4.1 (reduced Energy-Momentum Method). Suppose thatp x = FL(£<g(a;)) 
is a relative equilibrium of the simple mechanical system (|2.3p wii/i momentum /i 
and velocity £ € Assume that G M is compact, G x = H and that a (H n G M )- 
invariant splitting = g Px p. Lei f be the orthogonal projection of £ onio p. 
T/ien i/ 

(i) dim Q — dim G + dim G^ > 7 and 

(ii) (d^V^i + corr^j_ (x))\ s is positive definite 

or 

(i) dim Q — dim G + dim G x = ; and 

(ii) (d x V^± + corr^j_ (x))| s is definite (positive or negative) 
then the relative equilibrium is G ^-stable. 

Proof. According to Theorem l2.f I the relative equilibrium is stable provided d^/i^±| y 
is definite. As both K^± and V^± have a critical point at p x ((Hi) in Proposition 
I4.ip . and since h^± = K^± + Vc± , then 

With the isomorphism k in (|4.2|) we can compute each of these Hessians in V s 
parameterized by elements in q M S © S*. 

Let us compute first d^i^^^ . Let A, Ai,A2 S q M , a, a\,a^ £ S and j3, j3i,02 £ 
S*. Note that using (|4.2p we have that the vertical variation Sp v corresponding to 
an element A e q M , i.e. 7(<5p) = k (A, 0, 0), is 

5P V = (~P, [(DI.(e x ) Q (x))(A)+ad>],-ip s [(DI.(.))(£ x ,A)]V 
Using Lemma I4T21 we have 

^ y - • = (-p t [(di • aq(x))(£ ± )] , 0) . 

Similarly, if a e S we have for I(Sp) = n (0, a, 0) 

V = (-±P t [(DI.a)(^)],«G(a)(^),.» s ), 

and then 

<5p v - (T xX « X • <^) V = (-P r [(DI • a)(^)] , 0) . 
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Finally, if (3 £ S* for I(Sp) = k (0, 0, 0) we have 5p v = (0, /?) and Sp H = (0, 0) so 

Sp v -(T xX ^ ■Sp H ) v = (0,f3). 

According to (ii) in Lemma |4. II we can now write, for A, Ai, A 2 £ q^, a, a\, a 2 £ S 
and 0,0i,p2 G S* 

d^iQx (k(Ai, 0, 0), «(A a , 0, 0)) = (P r [(DI • A 1Q (x))(^)], H^P.KDI • A 2Q (:z))(^)]> 

d^iQx («(A, 0, 0), k(0, a, 0)) = (P r [(DI • A Q (a;))(£ x )],io ^[(DI ■ a)(^)]) 

d^iQx («(0, ai, 0), k(0, a 2 , 0)) = (P r [(DI • a 1 )(e ± )],I^ 1 P r [(DI • a 2 )(^)}) 

d2 x ^( K (0,0,/3!), K (0,0,/3 2 )) = «/3i,/3i»s- 

d2 x ^( K (0,0,/3), K (A,0,0)) = 

d 2 Px K^(K(0,0,/3), K (0,a,0)) = 0. 

Where <C •, • ^s* is the inner product in S* induced from -C •, • via the 
Riemannian Legendre map ()2.2p . We compute now the remaining contribution, the 
Hessian of the augmented potential energy. Using Lemma 14.11 it is immediate to 
obtain 

d% x Vgx («(Ai, ai, 0), k(X2, «2, 0)) = dj^x {\iq(x) + ai, A 2Q (a;) + a 2 ) 
d^^x(/s(A,a,0),«(0,0,/3)) =0 
d^^x(K(0,0,A),«(0,0, J 8 2 )) = 

for every A, AiA 2 £ q M , a, ai, a 2 £ S, /3, /3i, /? 2 € S*. 

Therefore, dp x /i ? x| y block-diagonalizes in the two blocks K(q M © S © {0}) and 
k({0} © {0} © S*) as 

dp^h^y (k(Ai,cji,0), k(A 2 , a 2 , 0)) = (d|Vgx + corr ? x (a)) (Aiq(x) + ai, A 2 q(x) + a 2 ) 
d2 x ^x| y S («:(A,a,0), K (0,0,/3)) =0 
dJ.VlviK ' '^ 1 )'^ ' '^)) = «/?i,A»s- 

for every A, AiA 2 £ q M , a, ai, a 2 € S, /?, 0i,02 £ S*. From the above expression and 
Definition 14. II it follows that the bilinear form d^ft^x^ is equivalent to the pair of 
bilinear forms (d^V^x + corr ^x (&))] and <C •, • 3>s* ■ 

Let us examine now when d^/i^x^ is definite. Since the block <C •, • ^s* is 
positive-definite or trivial, there are two possible scenarios, (i) dim S = and (ii) 
dim S > 0. Note that dim S = if and only if the dimension of the orbit G ■ x equals 
dim Q. Recall also that dim G ■ x — dim G — dim G x . In this case dp x h^±\ v consists 
only in the block (d^V^x + corr^x(a;))| s and the relative equilibrium is G M -stable 
provided this block is definite (positive or negative). 

For the other possibility, if dim S / then dim Q > dim G ■ x and the block 
■C •, • ^5>s' is positive definite since the metric on Q is Riemannian, so dp x /ijx| y 
is definite if and only if the block (d^Vjrx + corr^x (x))| E is positive definite. This 
completes the proof of the theorem. □ 

Remarks. 

• It is easy to see that the above theorem particularizes in the regular case 
(G Px = G x = 0) to the main result of the reduced Energy-Momentum 
Method, (see page 35 in [15]). In the regular case, one has to test the sta- 
bility of the Hessian d x V^ restricted to • x) ± , where p. is the momentum 
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value of the relative equilibrium under study and is Smale's amended 
potential energy (see [2D]), denned as 

V^{x) = V{x)+ l -(^l- 1 {x)(n)). 

Obviously this function is not well defined at those points x € Q such that 
dirng^ > 0, since the locked inertia tensor is not invertible. Hence, it is 
not possible to define the Hessian of the amended potential in the singular 
setting. However, the following relation is to be noted (see [19] . equation 
(2.28)), which holds in the regular case at a relative equilibrium 

d 2 Mv lt v 2 ) = d 2 x Vt( Vl ,v 2 ) + ((DI • wxXO,!- 1 ^) [(DI • v 2 )(Z)}). 

This suggests that what we are testing in Theorem l4.1l is exactly the singular 
analogue of the Hessian of the amended potential, even when we cannot talk 
about the amended potential itself. 
• In the regular case, according to our definition (|3.14|) of t, £ = (g M • x) 1 - 
with respect to -C •, • S>. In the presence of isotropy for the base point x of 
our relative equilibria, it can be seen that £ is the orthogonal complement 
to g M • x within the space of admissible variations (see (|7.1[) and the proof of 
Proposition 16 . 1[) . Therefore, conditions (ii) in Theorem 14.11 are tested in a 
space S which is orthogonally complementary (with respect to the kinetic 
energy metric) to the drift orbit G„ ■ x. 

5. An example: The sleeping Lagrange top 

We will apply our singular version of the reduced Energy-Momentum Method 
to the study of the relative equilibrium known as the sleeping Lagrange top. This 
problem has been extensively studied in the literature of Classical Mechanics; see 
in particular [8] for a geometric perspective of the problem using the Hamiltonian 
and symplectic formalism. Here we show the advantages of the method stated in 
Theorem 14.11 when studying stability in simple mechanical systems. Indeed, taking 
into account the extra cotangent bundle structure of the problem actually leads to 
simpler calculations for obtaining stability results when compared with the non- 
adapted methods constructed for general symmetric Hamiltonian systems. 

The Lagrange top is a symmetric simple mechanical system defined on the con- 
figuration space Q — SO (3), in the usual representation by orthogonal 3x3 real 
matrices with determinant 1, and being the symmetry group G = T 2 . We use the 
right trivialization for TSO(3) given by the isomorphism TSO(3) -> so (3) x SO(3), 
i.e. 5g = £g with £ G 9, and identifying so(3) with M 3 under the inverse of the usual 
isomorphism given by R 3 3 u — > u G so (3) (as a 3 x 3 matrix algebra), defined by 
uv = u x v, Vv G M 3 . Analogous considerations hold to obtain the trivialization 
of the phase space for the problem: T*SO(3) — > M 3 x SO(3). Then we can write 
the Hamiltonian for this system as 

(5.1) H(tt, A) = -7r • E^tt + mglAe 3 ■ e 3 , 

where A S Q, ir S K 3 , and [it, A) S T£Q, E = diag(i,i,i 3 ) and Ea = A£A*. 
Finally, to, <?, I are physical constants of the model. 

The group T 2 is identified with S 1 x S 1 where, in our matricial representation, 
both copies of the circle group act by rotations around e3 . Its action on Q is given 
by (L,R) - A = LAR* and the action on the phase space is by cotangent lifts. We 
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do not actually need the (simple) expression for the cotangent-lifted action since 
our methods rely finally on the geometry of the action of G on Q once they were 
constructed in this spirit. Let (I, r) e I 2 = Lie (T 2 ), then the infinitesimal action of 
R 2 on Q is given in the right trivialization of TSO(3) by (I, r)q(A) = (le 3 — rAe 3 , A). 

To finish with these preliminaries, let us show the Riemannian structure of Q 
associated to the kinetic energy of the problem. If V\ = («i,A) and V 2 = (1*2, A) 
are two tangent vectors to Q at A, then 

(5.2) < Vi, V 2 > (A) = vi ■ E A v 2 

which is easily checked to be a T 2 -invariant symmetric contravariant bilinear tensor 
on Q. Note also from the expression (|5.ip of the Hamiltonian, that the potential 
energy is given by V(A) = mglAe 3 ■ e 3 . 

It is well known that the phase space point pi = FL((Z, J-)q(I)), where I is the 
identity 3x3 matrix, is a relative equilibrium of this system, known as the sleeping 
Lagrange top. Our aim is to study its stability. First of all, note that since A = 7, 
then (Z,r)<g(7) = (C e 3i -0> where the number £ = / — r uniquely determines the 
infinitesimal generator, and it is physically interpreted as the angular velocity of 
the rigid body modelled by this system. 

The first thing we need to compute is the locked inertia tensor. It easily follows 
from the previous expressions that 

< (Zi,ri) Q (A), (Z 2 , r 2 )g(A) >= (he 3 - riAe 3 ) ■ E A {l 2 e 3 - r 2 Ae 3 ) 



= (/i,ri) 

and then 

1(A) = 



e 3 • E A e 3 -e 3 • AEe 3 \ f h 
-e 3 • EA l e 3 e 3 • Ee 3 J \ r 2 



e 3 • ^Ae 3 -e 3 • A^e 3 
-e 3 • £A*e 3 e 3 • Ee 3 

Thus at the configuration A = I of our relative equilibrium we have 

h -h 



Hi) , 

' -13 «3 

and so the momentum value of the relative equilibrium is 

H = !(/)(/, r)=t 3 (C,-C). 
As the configuration point of our relative equilibrium is the identity element of 
SO (3), then its isotropy group is H — Gi — S 1 , regarded as the diagonal embedding 
of S 1 in the 2-torus. For the momentum isotropy, just by noting that the group 
is Abelian and thus the coadjoint representation is trivial, we obtain G^ = T 2 . 
Finally, using the characterization G Px — G x n G M we also get G Vl = G x = S 1 . 

The next step is to choose a G Px -invariant complement to g Pa . in g M = R 2 in order 
to implement the splitting (13 . 1 3[) and obtain the velocity (I, r) . We will obtain 
it by orthogonality with respect to a G Px -invariant inner product in g M . Since the 
adjoint action is also trivial any inner product is ^-invariant. We can choose the 
family of inner products 

' k 
1 



G = 



Since g Px is generated by (1,1), its orthogonal complement p with respect to G 
is generated by the normalized vector k = , 1 (1, —k), so for the infinitesimal 
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generator of our relative equilibrium, £ = (I, r), we have that £ is the orthogonal 
projection of £ onto p, i.e. ^ = G(£,k)k = C (j^,-^)- 

The augmented potential energy V^± (A) = V(A) — |l(A)(^- L , corresponding to 
this relative equilibrium is now written as 

c 2 

V^± (A) = mglAe 3 ■ e 3 - 2 - ^ 2 [e 3 • £ A e 3 + 2fci 3 e 3 • Ae 3 + k 2 i 3 ] . 

We now compute the derivative of the augmented potential. According to our right 
trivialization, a tangent vector SA S TaSO(3) can be written as SA = eA, where e 
is a vector in R 3 . Then, after some vector calculus manipulations we get 
(5.3) 

c 2 

dV ( ± (A) • SA = mgle 3 ■ (e x Ae 3 ) - + ^ 2 [ki 3 e 3 ■ (e x Ae 3 ) - e 3 ■ E A (e x e 3 )] 

which vanishes at A = I, since FL((Z, r)g(J)) is a relative equilibrium. The Hessian 
of Fjx at A = 7 is 

dfV e x (M 1; M 2 ) = ( Cl x e 3 ) • (e 2 x e 3 ) ( ( *jff$ C ' - m 5 z) 
-ITw( ei xe 3)'^(e2 xe 3 ). 

A straightforward computation shows that the correction term (|4.4p vanishes. 

As the group is Abelian, then g M = g and q M = (0,0), so S = S. Then by 
Theorem 14. II the relative equilibrium is G M -stable if djV^x| g is positive-definite. It 
is then necessary to obtain the linear slice for the toral action on SO (3) at A = I 
with respect to the Riemannian metric (|5.2|) . A tangent vector SA G T/SO(3) 
written as SA = el is orthogonal to g ■ I if and only if 

e ■ Ee 3 = 

and thus we have 

S = {SA e T/SO(3) : e € span{(l, 0, 0), (0, 1, 0)}}. 

To check the positive definiteness of djV^±| s is then equivalent to showing that the 
2x2 matrix 

' (ki 3 + i 3 - i)( 2 
(1 + fc) 2 

is positive-definite, where I2 is the identity matrix in R 2 ; so the condition for 
stability is satisfied if 

(5.5) C 2 > g + k)2mgl 

ki 3 + i 3 - l 

Remark. The above expression is exactly the one obtained in [llj . page 718, by 
using the (singular) Energy-Momentum Method for general Hamiltonian systems 
in arbitrary symplectic manifolds, i.e. Theorem 1 2. II in this paper. In that work, the 
same condition for the Testability of the sleeping Lagrange top is obtained after 
computing algebraically the eigenvalues of a 4 x 4 matrix (d 2 ^/^^ ), which was not 
put in block-diagonal form by applying their general stability criterion. With our 
method the same expression follows easily from the unique eigenvalue of the scalar 
matrix (|5.4|) . This shows the potential power of employing methods adapted to the 
cotangent bundle structure in the study of simple mechanical systems. In particular 



M) rv:w* -mgi)h 



STABILITY IN SIMPLE MECHANICAL SYSTEMS 



17 



in practical stability problems with higher dimensional configuration spaces the 
implementation of Theorem 12.11 could involve checking the definiteness of much 
more complicated matrices, forcing the use of numerical methods in some situations 
where the reduced Energy-Momentum Method developed in this paper (Theorem 
14.11 and also Corollarv l6.2p could offer simpler or even exact results. 

For the sake of completeness we will sharpen the stability condition (|5.5| , follo- 
wing |llj . Since in the stability condition k appears, which is related to G, the 
sharpest (or optimal) stability condition will be the lowest value of C 2 among all 
possible k. This is an straightforward optimization problem in elementary calculus, 
so by differentiating the expression 

(i+^W 

kh + h - i ' 

we find that the minimum value is reached at k = 2 \~' 3 . So the sharpest stability 
condition yields the well-known lower bound for the angular velocity 

2 Amgli 
s ^ ■? ■ 

6. The singular Arnold form and Block-Diagonalization 

In the previous section we used the realization of the symplectic normal space V s 
at a relative equilibrium given in (|3.16[) and which has been shown to be isomorphic 
to q M © S © S* by the isomorphism (14. 2| . This was helpful to simplify the study of 
the definiteness of d^/i^^ , obtaining a block-diagonal structure with one block 
being positive-definite or trivial, reducing the problem to study the definiteness of 
(d^V^i + corr^i (x))| s . In this section we pursue the study of the symplectic nor- 
mal space V s by obtaining new block-diagonal forms for d^/i^l^ and also for the 
symplectic matrix fi at a relative equilibrium. Recall that in the /-representation, 
the symplectic matrix u>(p x ) of T Pn .(T*Q) has the following form (see [13] and [T7]). 
If (771,01; Pi, ax), (772,02; vi, oa) eI(T Px (T*Q)) = (rffiS)ffi (t* © S*) then 

(6.1) fi ((771,01; i>i,ai), (r7 2 ,a 2 ; f 2 ,a 2 )) = (^2,Vi) + ("2,011) - (^1,%) - (ai,o 2 ). 

In the following we will use a parametrization of V s different from (|4.2|) . This, when 
available, will present the extra advantage of putting d^/i^i| y in block-diagonal 
form consisting of three blocks instead of two, as in Theorem 14. 11 hence simplifying 
the stability analysis further. We will start by defining a singular analogue of the 
Arnold form, (see [1] and [9] to see how the Arnold form arises in the study of 
the stability of regular relative equilibria of simple mechanical systems). Hereafter, 
we assume that we will be working under the same conditions and hypotheses as 
in the previous section. In particular, we will always have p x = ¥L(£q(x)) as a 
relative equilibrium of the simple mechanical system (12. 3|) with momentum fi and 
base point isotropy G x = H . 

Definition 6.1. The singular Arnold form at a relative equilibrium with base point 
x and momentum fi is the bilinear form on q M , Ar : q p x q M — > K defined by 

Ar(Ai,A 2 ) = (ad^/j,A(x,/i)(A 2 )), 

where the map A(x,n) : q M — > r is defined by 

A(x, M )(A) = ig 1 (ad» +P t . [ad A (jfrV) 
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Remark. To the best of our knowledge, the first time a singular analogue of the 
Arnold form appears in the literature is in [6], equation (3.56), under the name 
of "generalized Arnold form" . In that work, this object is defined in the context 
of general Lagrangian systems and the Lagrangian Block-Diagonalization method. 
Also see Section [7] for a comparison of other results in [6] . 

We will now define another two spaces which will be useful for the obtention of 
block-diagonal expressions. The motivation for the introduction of these spaces is, 
following the ideas in |19) . that provided a non-degeneracy condition is satisfied, 
dp^/i^l^ will further block-diagonalized, refining the conditions of Theorem 14. II 

Recalling the identification k : q M S S* — > V s , let us define the following 
subspace of q M © S: 

(6.2) w mt = {(\ b ,b)eq»®S : (DI • (X b Q (x) + b))(^) G p*} . 

Also, using the map (A, b) i— ► Xq{x) + b which maps isomorphically q^ © S onto 
E C T X Q, wc define the following subspaces of E: 

(6.3) S rig = {Xq(x) : A S q M } , and 

(6.4) Stat = {\ b Q (x) + b : (A b , b) e Wint} • 
Note that we have the following obvious identifications: 

Srig ^ ^ 
Stat ~ Wtat • 

Following pi5] we can give the following interpretation for these two spaces: S is 
seen as the space of all admissible variations orthogonal to the infinitesimal drift 
directions g M • x. This space has a contribution E r i g which corresponds to variations 
in S which generate group motions, i.e. regarding our systems as a "rigid body" 
without internal structure. On the contrary, the subspace Ei n t corresponds to all 
the variations of our system that are purely internal, i.e. variations in "shape" , not 
coming from the symmetry group. 

Proposition 6.1. // the Arnold form is non- degenerate then E = E r ; g © E; nt . 
For the proof of this proposition, we will need the following lemma 

Lemma 6.1. For every A 6 q' 1 and v € S 

(i) P„ [(DI-Aq(x))(^)] =0 

(ii) P„ [(DI.«)(^)] =0. 

Proof. For (i), using (|3.3p . for any ( 6 f) one has 

(m-\ Q (x))(e,() = -I( 2; )(ad A C ± ,0-I(e ± ,ad A C) 
= -(ad AM ,C)=0, 

where the second equality follows since kerl(ir) = t) and I(x)(^ ± ) = /i. 

For (ii), making A = C € f) in item (iii) of Theorem l3.11 we have (DI • w)(^ ± , C) — 

for every £ e h. □ 

Proof, (of the proposition) It is clear that E r ; g C E and Ei nt C E so we have to 
prove S rig nStat = and S rig + E int = E. Let ^ A E q^. Then \q(x) S E rig nE int 
if and only if (DI(x) • Xq(x))(^ ± , e) = for every e e t + h. By (i) in Lemma loTTl 
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this holds if and only if the same condition is satisfied for every e £ i. Using (13. 3|) . 
this is equivalent to 

= I(o;)(ad A e ± ,e)+I(a ; )(e x J ad A e) = (I(a ; )(ad A (io 1 M))+ad>,e) 

= <i (P t . [ad A (to V)]) +ad A M,e) = (h (Pr- {(adA^ 1 m)] + I^(ad A /x)) , e) 

for every e£t, regarding that ad^/i G t if A G q M . Since lo is an isomorphism, this 
condition is the same as A(x,/i)(A) £ p C M , but then the Arnold form would be 
degenerate, which is a contradiction. 

To prove that E r j g + Ej nt — S let us note the following: if we call 

V={5q£T x Q : P^ [(DI-5g)(£ X )] = 0} , 

then by the definitions of t (|3 . 14[) . q M (|3 . 1 5|) . and by Lemma RTT1 we have that £ is 
precisely the orthogonal complement to g M ■ x in T>. The rest of the proof is then a 
consequence of Proposition 3.7 in [BJ. □ 

Corollary 6.1. If the Arnold form is non-degenerate then the map k : q^ 1 © Ej nt © 
S* — > V, defined as 

(6.5) «(A, (A&(s) + a), 7 ) = «(A + A a ,a, 7 ) 

/or every A € q M , (A a ,a) G Wi n t and 7 G S* is a G Vm -equivariant isomorphism. 

We will assume now that the Arnold form is non-degenerate and then we will 
study the symplectic matrix fl of V s and dz, x h^±. 

Proposition 6.2 (Block-Diagonalization forms) . If the Arnold form is non- degenerate, 
under the isomorphism k of Corollarv \6.1\ we have the following expressions for the 
symplectic matrix fl and dp x h^±\ v : 



(6.6) 



and 



(6-7) 2 1 

where the entries of f2 are 

S(Ai,A 2 ) = -(/i,ad Al A 2 ) 

*(A,(A^)+&)) = (/i,ad A A b ) 

^(( A Q^)+ a ) ) ( A Q( a; )+ & )) = -(A*,ad A£ ,A b )+«C( )(e ± ),a»s 

-«C(a)(£ x ),6» s ■ 

Proof. The form for fl follows trivially from (|6.ip and the definition of A; from 
(|6.5p . For dp^/i^^l^ , and recalling its block-diagonal form showed in the proof of 
Theorem 14. 11 the only two things that we must check are 

(i) (d£V£-L +corr 4 _L(x))(Ai Q (x),A 2 (5(a;)) =Ar(Ai,A 2 ), and 

(ii) (d 2 x V t ± +corr^(x))(\ Q (x),5q) = 0, 
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for A, Ai, A 2 e q M and Sq € Ei nt . 

To prove (i) recall that the potential energy V is G-invariant and then 

(d 2 x V^+corr^(x))(X 1Q (x),X 2Q (x)) = (-±d 2 x (l(x)(^, ^)) + corr^(^ (X 1Q (x), X 2Q (x)). 

We compute now both terms in the right hand side of the above expression. For 
the first one, 

-ldl(l(x)(e,e))(Xi Q (x),X 2Q (x)) = -±A 2Q ((DI.A 1Q (x))(£\n)M 

= A 2Q (I(ad Al £\^))(z) 
= (DI.A 2Q (x))(ad Ai e ± ,^) 
= -I^Xad^ad^),^) 
-I(:r)(ad Al ^,ad A2 ^). 

Now for the second, 

con^(x)(\ 1Q (x),\ 2Q (x)) = (P r [(DI-A 1Q (x))(^)] ,i - 1 P t [(DI-A 2Q (x))(^)]) 

= (I(.x)(ad Ai e ± ) +ad AlA1 ,i - 1 [l(x)(ad A2 C ± ) +ad A ^]) 
= I(.T)(ad A2 ^ ± ,ad Al ^ ± ) + (ad Al /j, ad Aa ^) 
+ (ad A2 /j,ad Al £ ± ) + (ad^^io^ad^/i)), 

where we have used that I(x)(?7i, 7y 2 ) = (77^ , 772) V 771 , 772 € 0, since kerl(x) = \). 
Also, (ad A /i,77) = (ad A /i,?7 r ) Vr; 6 9, since A S q^, which means that ad A /i 6 f)°. 
Finally, noticing that 

(ad A2 /a, ad Al ^) = I(x)(£\ ad A2 (ad Al ^)), 
and putting both contributions together, we obtain 
{dlV £ ± +cott c ±(x))(X 1q (x),X2q(x)) = (ad Al /j,ad A2 ^ +I ( 7 1 (ad A2 /j)) 

= (ad Ai /i,P r [ad A2 (to V)] +i 1 (ad A2 /j)) 
= Ar(Ai,A 2 ), 

since ad Ai /j £ \)° and ^ = l^ 1 /1. 

Again for (ii) because the potential energy V is G-invariant we have 

(d 2 x V ( ± +corr^(x)) (X Q (x), X a Q (x) + a) = 
= (-§d» +corr^(a:)) (A Q (x), A Q (z) + a). 

We will start by computing the contribution of the correction term: 
corr C x (») (A, A Q (*) + a) = (P r [(DI • A Q )] , jfr *P t [(DI ■ (A Q (*) + 0)) )] > 

= -<P t [lOrXadA^)] +ad>,i- 1 P r [(DI.(A«( 2 ;)+ a ))(^)]) 
= - <!(.*) (ad A ^) + ad A /i,I^(DI • (X a Q (x) + a)(^)), 

since kearl(z) = f) and by Lemma [6J] Pf, [(DI ■ (A|(i) + a)(^)] = 0. For the first 
term we have 

-~d^ , e)) (A, A Q (x) + a) = (DI • (AfcOr) + a))(ad A £ x , 

and so we finally obtain 
(d^x +oorr f x(a)) (Aq^), A Q (z) + a) = -(ad^I^DI ■ (A Q (.x) + a)(£ x )). 
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This expression is zero since by construction (DI • (Ag(x) + a))(C ± ) annihilates t if 
(A a , a) G ui; nt and I ( J" 1 (ad^/i) G t for every A 6 q M . To see this, note that the image 
of Ig 1 is in t and hence I(x)(Iq 1 (ad^/z), Q = (ad^/i, () — — (ad^, A) = 0, for every 
C G g M , in particular if £ G p. □ 

An inspection of the form of d^/i^i| y in (|6.7p together with Theorem 14.11 leads 
to the following sharper result concerning the G^-stability of p x . 

Corollary 6.2 (Block-diagonalization and stability). In the hypothesis of Theorem 
\4-l\ and assuming that the Arnold form in non-degenerate, if: 

(i) dim Q - dim G + dim G x > 

(ii) (d x V^± + corr^i (x))| s ^ is positive definite 

(iii) the Arnold form is positive definite, 

or 

(i) dim Q - dim G + dim G x = 

(ii) the Arnold form is definite (positive or negative), 
then the relative equilibrium is G^- stable. 

Remark. In [19], Theorem 2.7, equivalent results to our Proposition l6 . 21 and Corol- 
lary [621 are obtained as a consequence of their reduced Energy-Momentum Method, 
for the particular case G x = {e} (regular relative equilibria). 

7. Some remarks on the stability results 
The residual symmetry sub-blocking. 

In this subsection we will use the fact that the symplectic normal space V s supports 
a linear representation of G Px in order to improve our block-diagonalization results. 
We start with the following lemma (see [T7] for a proof): 

Lemma 7.1. Let h G H = G x , £ G r and v, w G T X Q, then 

(i) Adl-r [(DI • v)(Q] - (DI • (h ■ v))(Ad h 

(ii) < C(h ■ v)(Ad h £),h ■ w > s =< C(v)(0, w >s- 

As an immediate consequence of this and just by regarding their definitions (14. 3|) . 
(|3.15|) . (|6.4|) . the spaces Sjq^ and E; nt are G Px -invariant. We will use a tool from 
representation theory known as the isotypic decomposition of a linear space acted 
linearly upon a compact Lie group to take advantage of the residual symmetry 
group G Pt in order to further block-diagonalize d^/i^x]^ . For this, we need the 
following definition, which is taken from [3]. 

Definition 7.1. Let K be a compact Lie group acting linearly on a (real and 
finite dimensional) linear space N . The isotypic decomposition of N is the unique 
decomposition 

N = Ni®...®N r , 

where each Ni is the direct sum of all K -isomorphic irreducible subspaces of N, and 
it is called an isotypic component of N. 

The isotypic decomposition of a linear space satisfies the following remarkable 
property (see [3]): If B is a iC-invariant bilinear form on (that is, B(g-vi, g-v%) = 
B(v\,V2) for every g G K, v-y , t>2 G N) then B\ N x n = for every pair of isotypic 
components of TV with Ni ^ Nj . Therefore the expression of B block-diagonalizes 
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with respect to the isotypic decomposition of N. We will apply this property to 
the bilinear form given by d x V^± + corr^i(x). 

Lemma 7.2. Let p x be a relative equilibrium of the simple mechanical system (|2.3p 
with velocity £. Then the bilinear form d x V^± -I- corr ^ (a;) is G Px -invariant. 

Proof. We have to prove separately the invariance of each term. For d^V^j_ the 
result follows if we prove that V^± is G Px -invariant. Since V is G-invariant, and 
so G Px -invariant, we only need to prove that the function §I(')(s j£ ) is G Px - 
invariant. Recall that at a relative equilibrium, besides the characterization G Px = 
G M n G x , one also has G Px = {h £ G x : Adh^ = £ }, which follows trivially from 
the property 

g £ G Px «■ g ■ FL(^(z)) = FL(#(x)), 

since £ G p C t and by definition our relative equilibrium can be written as 
p x = FL(£q(.t)) = FL(£q(x)). Recall also the invariance property of the locked 
inertia tensor (|3.2p . Then, for every x' £ Q and h £ G Px , one has 

h(h-x')(e,^) = h(x')(Ad h -^,Ad h -^) = h(x')(^,e)- 

For the correction term, recall from Lemma 17. 1[ that if h £ G Px and 5q £ T X Q, 
then 

(DI • (h ■ Sq))^) = Ad^-i ((DI ■ 5q)(^)) , 

since Ad/i£ = Note also that to : r — > r* is a G x -equivariant isomorphism, 
that is Iq o Adh — Ad^-i o I . Then, given Sqi, 8q2 £ T X Q and h £ G Px , we have 

COTT^(x){h ■ Sq u h ■ Sq 2 ) = ((DI • (h ■ SqxM^X 1 ((DI • (h • MX^))) 

= (Ad*-, ((DI • S qi )(e)) X 1 (Ad^-x ((DI • Sq 2 )(e)))) 
= (Ad^-x ((DI • ^)(^)) , Ad h (V ((DI • 5q 2 )(^)))) 

= ((Di-<5 (Zl )(^),i ( 7 1 (( DI -^2)(C ± ))) 

= corr C i(a;)((5(ji,(5g2) 

□ 

Thus, in the hypothesis of Theorem 14.11 we can use the fact that £ is a G Px - 
invariant subspace of T X Q, and then testing the definiteness of (d^.V^± +corr ^± (x))\ s 
is equivalent to testing the definiteness of every restriction (d x V^± + corr^ (x))| E 
where 

£ = o . . . O S r 

is the isotypic decomposition of S. Analogously, if the Arnold form is not de- 
generate, to test the conditions for stability of Corollary 16.21 is equivalent to test 
definiteness of Ar q f and (d x V^± + corr^i (x))| s for each of the isotypic compo- 
nents of q M and £j nt , respectively. 
Nature of the stability results. 

In Theorems 12.11 and 14.11 we have imposed the compactness of the momentum 
isotropy subgroup G M . However, this compactness condition can be weakened by 
assuming that y, is split, and this is how the original theorems are stated in [5] and 
[llj . A momentum value [i is called split if there exists a G^-invariant complement 
to g M in g. Obviously this is the case if G M is compact, since in this case one can 
define this complement to be the orthogonal complement to with respect to some 
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G^-invariant inner product on g. Likewise, if G itself is compact or Abelian, then 
every momentum value is automatically split. 

In the most general situation, if the relative equilibrium under study has not 
a split momentum value, Theorem 14.11 and Corollary 16.21 are still applicable, but 
in that case one does not obtain conditions for G^-stability, only for the weaker 
notion of leafwise stability. A relative equilibrium is called leafwise stable if it is 
G^-stable for the restriction of the Hamiltonian flow to J _1 (/i). The reason for this 
nomenclature is that in the free case, a relative equilibrium z with momentum fi 
for the symmetric Hamiltonian system (V, u), G, J, h) is leafwise stable if the point 
[z] £ V/G in the orbit space is a Lyapunov stable equilibrium for the reduced 
Hamiltonian system on the symplectic leave J _1 (/i)/G A1 C V/G, rather than being 
stable in the full Poisson quotient V/G. The results of [3 [TTJ [15] guarantee that 
if /i is split, then leafwise stability of z implies Lyapunov stability of [z] in V/G, 
and, hence, so do Theorem 14.11 and Corollary 16.21 See [TB] for a more detailed 
explanation of these concepts. 

The reduced Energy-Momentum Method and Lagrangian Block- 
Diagonalization. 

In [6] a method was constructed for testing the stability of relative equilibria of sym- 
metric Lagrangian systems. In that work, the techniques of the reduced Energy- 
Momentum Method of [TH] are translated to systems defined on the tangent bundle 
TQ of the configuration space and developed for general Lagrangian systems in- 
variant under a possibly non free, tangent lifted action. We briefly explain the 
relationship of the results of [6] applied to simple mechanical systems with our 
work. See j6[ [7] for more details on the Lagrangian Block-Diagonalization method. 

Let L £ G G (TQ) be a function on the tangent bundle of Q invariant under 
the tangent lift of a proper action of the Lie group G on Q. This function is 
called a Lagrangian. There is a well-known procedure to obtain a bundle map 
FL : TQ — > T*Q constructed from L (no Riemannian structure is in principle 
available on Q). In the case FL is a diffeomorphism, the Lagrangian is called 
hyper-regular, and one can pull-back the canonical symplectic form from T*Q to 
TQ and define a Hamiltonian system on TQ (see [6] for details). Given an element 
£ € 0, the locked Lagrangian L% £ C°°{Q) is defined as 

h(x) =L(£ Q (x)). 

Also, the locked momentum map is defined as the map 1^ : Q — > g* that satisfies 

for every r\ £ g. We will also need the definition of the space of admissible configu- 
ration variations at a point x £ Q, which is 

(7.1) V = {Sq £ T X Q : D/ 5 • Sq £ °}, 

for a fixed element £ £ g. Finally, given x £ Q and £ £ q the linearized momentum 
map I x : q — ► g* is defined as 



UV) = 4 
y " dt 



t = W*) 



for every £, r\ £ q. We will denote its generalized inverse by I x 1 : range 
0/kerI x . 
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The Lagrangian Block-Diagonalization method gives sufficient conditions for for- 
mal stability of relative equilibria in Lagrangian systems. Under some assump- 
tions, formal stability implies G^-stability. Here we shall not be concerned with 
those differences, in order to provide a comparison between the Lagrangian Block- 
Diagonalization and our singular reduced Energy-Momentum Method. Further- 
more, one needs a technical condition relating 1^ and I x at a relative equilibrium in 
order to be able to apply the method (see (3.15) in [6]). However, for the particular 
case of Lagrangian systems defining simple mechanical systems this is automati- 
cally satisfied and thus it is not necessary for our comparison objective. So we 
assume that the Lagrangian Block Diagonalization method is applicable in order 
to simplify the exposition. 

The next proposition collects the results of [6] concerning relative equilibria and 
their G^-stability. Let g be the symmetric bilinear form on T X Q induced by the 
hyper-regular Lagrangian L and defined by 

d 

■ ods 

By the hyper-regularity hypothesis of L, g(-, •) is non-degenerate. 

Proposition 7.1 (Lagrangian Block-Diagonalization, [6]). Let L E C G (TQ) be a 
hyper-regular Lagrangian invariant under the tangent lifted action of the Lie group 
G on Q. A point v x S TQ is a relative equilibrium for the Lagrangian system 
defined on TQ by L if there is an element £ G Q such that v x = £q(x) and x is a 
critical point of L^ . If the bilinear form 



i \ d 
g(v x ,w x ) = — 



L(tv x + sw x ). 



B= ( -d*L € + <D/ 4 ■(■), i*" 1 (DJ f •(•))) 



V 



is positive (negative) semi-definite with kernel g^ ■ x and gi( fl . x )-L is positive (nega- 
tive) definite then the relative equilibrium is G ^-stable. 

We now study Lagrangians defining simple mechanical systems. It is a standard 
fact that if we are given a G-invariant Riemannian metric C > on Q and a G- 
invariant function V € C G (Q), then the Lagrangian formulation of the associated 
simple mechanical system (|2.3p is 

(7.2) L(v x ) = i || v x || 2 -V(x). 

For L of the form (|7.2j) it is straightforward to compute 

h = Kx) 
g =<•,•>. 

Therefore, relative equilibria of the simple mechanical system (|7.2I) are defined by 
a velocity ( £ g and a critical point x € Q of recovering the well- know result 
for critical points of the augmented Hamiltonian in simple mechanical systems. 

To see that the stability conditions of Proposition 17.11 are then equivalent to 
those given in Theorem 14. 11 one only needs to prove that at a relative equilibrium, 
the space S c T X Q of Definition ^. H is indeed a complement to g M • x in T>, but this 
follows from Lemma l6.1l and it has been already used in the proof of Proposition [6Tl 
This shows that our stability result, Theorem 14. 1( is a Hamiltonian version of the 
Lagrangian Block-Diagonalization method applied to simple mechanical systems. 
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In the same way, one can see that the extra block-diagonalization construction 
carried out in subsection 3.3 of [6] is a consequence of the splitting of in rigid 
and internal subspaces of Section [5] in this paper. 
The normal form for the symplectic matrix. 

Besides the convenient form for d^/i^l^, given in Proposition ^. 21 which allowed us 
to refine Theorem 14. II and re-express it as Corollarv l6.21 it is important to note the 
particular expression for the symplectic matrix Q of the symplectic normal space 
V s identified with q M © £i nt © S* . In the free case (see [9] and [19]), the explicit form 
for the symplectic matrix, together with the one for d^/i^ly , allows the authors 
to obtain the linearization of the Hamiltonian vector field at a relative equilibrium. 
This is an important observation, since the study of this linearized vector field 
has applications in the spectral and linear stability (and instability) of the relative 
equilibrium under study, as well as for the identification of possible bifurcations 
from it. 

It seems that despite the generalization to the Lagrangian side and non-free 
actions in [6] of the stability results provided by the reduced Energy-Momentum 
Method of [19] , this feature has not been studied in detail for group actions with sin- 
gularities and the expression for the symplectic matrix obtained in Proposition 16. 21 
for relative equilibria with non-discrete configuration isotropics cannot be found 
in the literature. Here we prove that our expression (16. 6|) coincides in the regular 
(free) case with the one obtained in [15], equations (2.83) and (2.85). 

Indeed (and if we assume that the Arnold form is non-degenerate) given two 
elements (»7i, (Ag*(af) + a;), ft) € q M © lOint © S* ~ V a , i = 1,2, we will write 
5qi = (\ a *) Q (x) + <n and a(6qj = I^ 1 (J(FL(<%))) = A a \ Then we have the 
following: 

Proposition 7.2. The expression for the symplectic matrix Q of Proposition [fOl 
is equivalent to 

n((rn,6qi, (3i), (r) 2> 5q 2) fo)) = (^-[vuVs] ~ [Vi) a ( s Q2)] + [%> "(<%)]) 

+(ft>,<%) - (Pi,Sq 2 ) - dx e {x)(Sq 1 ,Sq 2 ), 

with defined in (|4.1[) . This coincides in the regular case with equation (2.83) 
in QJ]. 

Proof. Recall from Proposition 16.21 that we have 

n(( m ,6q ll p 1 ),(r l2 ,6q 2l p 2 )) = {fi,—[rji, r) 2 ] — [r)i, A" 3 ] + [772, A" 1 ]) 

+{P 2 ,a 1 )-(0 1 ,a 2 )-( f x, [A ai , A" 2 ]) 
+ < C(a 2 )^ ± ),a 1 >s - < C(ai)(£ x ),a 2 >s • 

Now note from the definition of a, that [r]i, a(Sqj)] = [i]i,\ aj ]. Also, since ft € 
S* = (g • x)° , i = 1,2, one has that (ft, 5^-) = (ft,a,-). Therefore, the proposition 
will be proved if we show that 

dx^(x)(S qi ,dq 2 ) = ( M , [A ai , A a2 ])— « C(a 2 )(^), ai » s + « C( ai )(^),a 2 » s . 

To see this, we choose local extensions X — (A Qi )q + ai and Y = (A° 2 )q +a 2 of 
Sqi and 5q 2 near x and then we use the formula for the exterior derivative 

d X ^ (X, Y) = X ( X «" (Y)) - Y (x 6± (X)) - x^ ([X, Y}). 
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We will need the following lemma, proved in [17] , which shows additional properties 
of the family of local vector fields introduced in Section 3. 

Lemma 7.3. The local vector fields defined by (|3.11[) satisfy 

(i) [v a ,Vb] (x) = for every v a , v b e S 

(ii) [Aq,tJ] (x) = for every A 6 r, v € S 

(iii) There is a small enough open neighbourhood O 3 x such that for every 
x'eG-x(~]0,v<ES and £ € g, < v(x'), £q(x') >= 0. 

(iv) Let x 1 — o~([g, s]) with a defined in (|3 . 10[) and A G r, then 

Xq ^ = ^ t = CT ([( ex Pe iA )5' s ])- 

Recall the definition \^ = ^L(^q) and let us compute separately the group 

orbit and slice contributions for dx^ (x)(X,Y) using the usual properties of the 
Levi-Civita connection. 

dx $ (x){a 1 ,a 2 ) =a 1 (< £^,a 2 >) (x) - a 2 (< £q-,3i >) (x) 
=« V Sl ^(x),a 2 > + < ^q(x),V Si S 2 (x) > 
- <C Va 2 ^(a;),ai > - < £q- (x), Va 2 ai(x) > 
=< C(ai)(^),a 2 >s - < C(a 2 )(^),ai >s, 

where we have used the definition of C (|3.12[) , together with item (v) in Theorem 
13.11 and the fact that 

<^g(x),Va 1 a 2 (x) > - <£Q(x),Va 2 ai(a;) > 

=< £g(x),T(ai,a 2 )(x) > + < ^(x), [ai,a 2 ](x)] >= 

by item (i) of Lemma 17.31 as well as noting that the Levi-Civita connection has zero 
torsion (T(X, Y) = V X Y - VyA - [X, Y] = 0, for every pair of vector fields X, Y). 
In the same way we can compute 

d x ^(x)((X^) Q (x), aj ) = (A q <)q (« #,a, ») (x)-aj (« ^,(A°')q ») (*) 

([(A°*) ,^]) (x) = -(pi- aj )(e,X ai ), 
since by (iii) and (ii) in Lemma 17.31 the first and last contributions vanish. Finally 

dx« i (x)((A« 1 ) Q (x),(A^) Q (x)) = 

(A Q1 )Q (« (A" 2 )q ») (x) - (A° 2 )q (« (A Qi )q ») (x) 
+ « # (x) , ( [A<» , A Q2 ] ) g (x) »= (DI • (A" ) Q (x)) (^ , A° 2 ) 
-(DI • (A Q2 ) Q (x))(e ± , A^) + (a, [A ai , A° 2 ]). 

Putting together all the contributions we obtain 

d x t ± (x){5q 1 ,6q 2 ) = (/i, [A Q1 , A Q2 ])+ « C(ai)(£ x ),a 2 » s - « C(a 2 )(^),ai »s 
+ (DI • Sqi)^, A a2 ) - (DI • 6&)(£ x , A ai )- 

The last two terms of the above expression vanish since for i = 1,2, A Qi e f C t 
and (A a \ a,) e w int , and by its definition ([621) p t [( DI ' <*9i)(£ X )] = 0. This finishes 
the proof. □ 
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